Abstract. We have performed new determinations of sul- 
phur and silicon abundances for a sample of 26 disk stars 
based on high-resolution, high signal-to-noise spectra. The 
results indicate a solar [S/Fe] for [Fe/H] > —0.3, below 
which [S/Fe] increases to ~0.25 dex at [Fe/H] = -1.0. We 
find that there is a good correlation between [S/H] and 
[Si/H], indicating the same nucleosynthetic origin of the 
two elements. It seems that the ratio of sulphur to silicon 
does not depend on metallicity for [Fe/H] > —1.0. The 
implications of these results on models for the nucleosyn- 
^vq [ thesis of a-capture elements and the chemical evolution of 
■ the Galaxy are discussed. 
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1. Introduction 

The trend of [O/Fe] vs. [Fe/H] has recently been paid 
much attention because works by different authors yield 
very dissimilar results. Specifically, fsraelian et al. ( |1998| , 



2001) and Boesgaard (199£) derived oxygen abundances 
from OH UV lines and the infrared O I triplet and found 
a linear rise of [O/Fe] to 1.0 dex when the metallicity 
decreases to [Fe/H] = —3.0 dex. Some doubts have, how- 
ever, recently been cast on the correctness of these re- 
sults from calculations based on the new generation of 3D 
hydrodynamical model atmospheres (Asplund & Garcia 
Perez 2001 ). Other studies (e.g. Nissen et al. 1992 , Nissen 
et al. 2001 ) based on the forbidden line at A 630.03 nm 
suggest, on the other hand, a more or less flat trend of 
[O/Fe] in the metallicity range of -2.0 < [Fe/H] < -1.0. 
Concerning other well-studied a-capture elements, e.g. 
Mg, Si and Ca, earlier observed trends of [a/Fe] versus 
[Fe/H] (e.g. Ryan et al. 1996D seemed to be in good agree- 
ment with a plateau of [a/Fe] ~ 0.4 dex. In particular, 
Fuhrmann (1998) found a clearly flat [Mg/Fe] ratio for 
[Fc/H] < —0.8 based on a LTE abundance analysis of Mg 
and Fe. However, the recent work of Stephens & Boes- 
gaard ( |2002 ) based on Keck spectra of more than 50 halo 
and disk stars suggests a quasi-linear trend of increasing 
[a/Fe] with decreasing [Fe/H] (see their Fig. 20), but we 
note that their average [a/Fe] is only about ~ 0.4 dex at 
[Fe/H] = -3.0. In addition, Idiart & Thevenin ( |2000D have 
found a complicated structure in [Mg/Fe] and [Ca/Fe] 
vs [Fe/H] based on a non-LTE re-analysis of equivalent 
widths from various sources. Although the scatter is sig- 
nificant for [Fe/H] < —1.0, their results are consistent with 
plateau-like mean trends at ~ +0.2 — 0.3 for both [Mg/Fe] 
and [Ca/Fe] for halo stars. Since different trends of [a/Fe] 
vs. [Fe/H] correspond to different nucleosynthesis histo- 
ries of the a elements, a clarification of the divergency on 
the trend of [O/Fe] vs. [Fe/H] is very important for many 



astrophysical fields including stellar evolution theory, nu- 
cleosynthesis theory, the formation and chemical evolution 
of the Galaxy, and the age of the Galaxy and the Universe. 

Another a element, sulphur, is highly interesting as an 
independent tracer of the nucleosynthesis of a elements 
and a well-established trend of [S/Fe] vs. [Fe/H] could be 
helpful in clarifing the discrepancy of the [O/Fe] ratio. Fur- 
thermore, sulphur is of particular interest in connection 
with studies of the chemical enrichment of damped Lya 
systems being one of the few elements that is not depleted 



onto dust (Centurion et al. 2000). Unfortunately, sulphur 
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has been ignored for many years probably due to the lack 
of measurable lines in the optical spectra. Presently, our 
knowledge of sulphur abund ances is limi ted to only a few 
studies: Clegg et al. (|l98l|) , Frangoi s (|l987t |l988|) and 
recent work s by I sraelian & Rebolo ( 2001 ) and Takada- 
Hidai et al. ( 2002 ). Again, the disagreement on the [S /Fe ] 
ratio at low metallicities appears: Frangois ( 1987 , 1988| ) 
suggested the [S/Fe] ratio to be constant for star s with 
-1.6 < [Fc/H] < -1.0, while Israelian & Rebolo fl200l| ) 
favor a linear rise of [S/Fe] ratio with decreasing [Fe/H] 
in the metallicity range of -3.0 < [Fe/H] < -0.6. With 
a smaller slope (—0.25 vs —0.46 in Israelian & Rebolo 
paper) the work of Takada-Hidai et al. ( 2002| ) support a 
linear trend of [S/Fe]. 

Sulphur and iron abundances in disk stars are usually 
better determined than in the case of halo stars. The sul- 
phur lines have a suitable strength for accurate abundance 
determination and Fe n lines, which have very small non- 
LTE effects, can be used for deriving the iron abundance. 
Therefore, the [S/Fe] vs. [Fe/H] trend can be more reli- 
ably established. With this advantage, the determination 
of sulphur abundances in disk stars with [Fe/H] > —1.0 
is important especially from the nucleosynthesis theory 
point of view. The comparison of the sulphur behaviour 
with those of other a elements in disk stars provides a key 
clue to the genuine origin of sulphur and further to the 
chemical evolution of the Galaxy. As the nearest even-Z 
neighbour of sulphur in the periodic table of elements, sil- 



Y.Q. Chen et al.: Sulphur abundances in disk stars: a correlation with silicon 



3 



icon is very suitable as a representative of the a elements 
for such a comparison. 

In the present work, we investigate the behaviour of 
sulphur in disk stars and the nucleosynthesis origin of this 
element. The investigation takes benefit of more accurate 
observations of sulphur lines, better determinations of stel- 
lar parameters, improved atmospheric models, and hence 
more reliable derivations of sulphur, silicon and iron abun- 
dances than in previous works (Clegg et al. 1981; Francois 



1987 , 1988 ), where the scatter of [S/Fe] for disk stars is 



so large that it is not clear if sulphur really behaves as 
an a-capture elements or not. In a following paper (Nis- 
sen et al. 2002), the study will be extented to halo stars 
such that the trend of [S/Fe] will be covered for the whole 
metallicity range of -3.0 < [Fe/H] < +0.2. 



2. Observations 

The stars were selected from the uvbyft photometric cat- 
alogues of Olsen (1983, 1993) according to the criteria 



of 5600 < T cff < 6600 K, 3.8 < log g < 4.5 and 
— 1.0 < [Fc/H] < +0.5 with approximately equal num- 
bers of stars in every metallicity interval of 0.1 dex. This 
selection aimed at obtaining a sample of stars with similar 
parameters in order to minimize the analysis uncertainties 
including the ways of determining stellar parameters and 
atmospheric modeling. 

Spectra of stars were obtained during three observ- 
ing runs: Oct. 17-22, 2000, Jan 8-10, 2001 and Aug. 11, 
2001 with the Coude Echelle Spectrograph attached to 
the 2.16m telescope at National Astronomical Observato- 
ries of Chinese Academy of Sciences (Xinglong, China). 
The detector was a Tek CCD (1024 x 1024 pixels with 
24 x 24 nm 2 each in size) . The red arm of the spectrograph 
with a 31.6 grooves/mm grating was used in combination 
with a prism as cross-disperser, providing a good separa- 
tion between the echelle orders. With a 0.5 mm slit (1.1 
arcsec), the resolving power was on the order of 37 000. 
The exposure time was chosen in order to obtain a signal- 
to-noise around 150-250 and the spectral coverage is 580- 
880 nm except for 3 stars with the wavelength range 550- 
830 nm. As examples, three sections of the spectrum for 
HD 168151 and HD 22879 are shown in Fig. and Fig. | 
with Fell, Si I and Si lines marked by vertical lines. Note 
that HD 22879 is the most metal-poor star in this sample, 
and the weak Si 869.39 line was not used in deriving the 
sulphur abundance because the noise level in the contin- 
uum at this region is comparable with the line strength. 

The data reduction, using MIDAS software, involved 
standard routines for order definition, background sub- 
traction, wavelength calibration, flat-field division, spec- 
trum extraction followed by radial velocity shift correc- 
tion and continuum normalization. The equivalent widths 
(EWs) were measured by direct integration or Gaussian 
fitting depending on which method gave the more reli- 
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Fig. 1. Three sections of the spectrum for HD 168151 
(T eff = 6530 K, log g = 4.12, [Fe/H] = -0.28) as observed 
with the 2.16m telescope at Xinglong Station. 



able value. The two Si lines at A869.4nm were separately 
measured using a combination of two Gaussians. 

Spectroscopic binaries and a few stars with a high- 
rotation of K-otsinz > 25kms _1 were excluded from the 
sample leaving finally 26 stars for abundance analysis. Ta- 
ble [l] presents the star sample as well as the stellar param- 
eters, kinematics, the equivalent widths of the S I lines and 
abundances. 



3. Analysis 

3.1. Stellar parameters 

The effective temperature was estimated from Stromgren 
photometric indices (Olscn 1983 , 1993 ) by using the 6 — y 
calibration of Alonso et al. (1996). We have calculated the 



E(b — y) for every program star using the calibrations of 
Olsen ( 1988 ) but found no significant reddening correc- 
tion. All stars have E(b — y) less than 0.010. 

K magnitudes for 4 stars are available in the 2MASS 
survey (Finlator et al. 200C| ). We derived the tempera- 
tureby using the V — K calibrations of Alonso et al. 
( |1996| ). The results are 5698 K for HD 22879, 6077 K 
for HD 28620, 6232 K for HD 77967 and 5882 K for 
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Fig. 2. Three sections of the spectrum for HD 22879 
(T off = 5788 K, log g = 4.29, [Fe/H] = -0.91) as observed 
with the 2.16m telescope at Xinglong Station. 



HD 216106. These values are slightly lower than those de- 
rived from the b — y indices (see Table [j]) with a mean 
deviation < T(V - K) - T(b - y) >= -66±27 K, but 
they are consistent within the estimated error of the tem- 
perature determination, i.e. ~100 K. For consistency, we 
adopted the value from the b — y indices for all stars. Fur- 
thermore, we checked that there is no significant trend of 
the derived iron abundance with excitation potential of 
low energy level of lines, supporting the adopted temper- 
atures to within ±100 K. 

As the gravity determines the degree of ionization for 
each element and it has a significant effect on Fe n and S I 
lines, a careful derivation of this parameter is important. 
With the available Hipparcos parallaxes, we can deter- 
mine accurate gravities from the relations g oc A4 / R 2 and 
L oc i? 2 T e ff 4 . The stars have accurate parallaxes with a 
relative error of < 0.05 except for 3 stars with ^ ~ 0.08. 
Using the derived temperature and absolute magnitude, 
the mass was estimated from the stellar position in the 
My — log T c ff diagram by interpolating in the evolution- 
ary tracks of VandenBerg et al. ( 200C ) and the bolometric 
correction was interpolated from the grids of Alonso et al. 

( HP - 

The metallicity, required in the input for temperature 
determination and abundance calculation, was first de- 
rived from the Stromgren rn\ index using the calibrations 



of Schuster & Nissen (1989). But the final metallicity from 
the abundance analysis of the Fe n lines was used to iterate 
the whole procedure to full consistency. The microturbu- 
lence, £ t , was obtained by requiring a zero slope of [Fe/H] 
vs. EW based on Fe II lines. 

The uncertainties of the stellar parameters are esti- 
mated to be around 100 K in T c g, 0.1 dex in log g, 0.1 
dex in [Fe/H] and 0.3kms _1 in £<. 

3.2. Abundance calculations 

The abundance analysis is based on a grid of flux constant, 
homogeneous, LTE model atmospheres , which were com- 
puted with a recent version of the MARCS code using 
updated continuous opacities including UV line blanket- 
ing of millions of absorption lines and many molecular 
lines (Asplund et al. |1997 ). These models are based on 
the assumption of [a/Fe] = 0.0. Differences are only 0.05 
dex in [Fe/H], -0.04 dex in [Si/Fe] and -0.02 in [S/Fe] by 
using a model calculated with [a/Fe] = 0.4 for HD 22879, 
which is the most metal-poor star in our sample. Given 
that enhancements of [S/Fe] and [Si/Fe] for our sample 
of stars are less than 0.3 dex, the enhancement of [a/Fe] 
or not in the model atmosphere has a negligible effect on 
the result. The abundance analysis program, EQWIDTH, 
was used to calculate the theoretical equivalent widths 
from the models, and elemental abundances were derived 
by requiring that the calculated equivalent widths should 
match the observed ones. When more than one line of the 
same species were measured for a star, the mean value is 
adopted by giving equal weight to each line. 

Loggf values for the Si, Si I and Fen lines are pre- 
sented in Table || together with solar equivalent widths 
measured from the Moon spectrum observed at Xinglong 
except for three lines with A > 860 nm. The equivalent 
widths for these three lines were measured from the so- 



lar flux spectrum of Kurucz et al. (|1984|) because the 
Moon spectrum observed at Xinglong does not cover the 
wavelength range of A > 820 nm. Solar abundances are 
calculated based on a solar model atmosphere computed 
with the same code as used for the stars. The (//-values 
from Biemont et al. (1991) (including the normalization of 
—0.05 dex to get agreement with lifetime measurements) 
are adopted for most Fen lines. With these values, we de- 
rived a solar Fe abundance close to the meteorite value 
loge Q (Fe) = 7.50 (Grevesse & Sauval |l99Sj ). We obtained 
differential \oggf values for the remaining three Fen lines 
by forcing them to give a consistent solar abundance. The 
logg/ values for four of the five S I lines are taken from Ku- 
rucz ( 2002 ) . As shown in Table |2|, the line-to-line scatter of 
the derived sulphur abundance is small for these four lines. 
From the A604.60nm line, we get a higher solar sulphur 
abundance. In the solar flux spectrum, this line seems to 
be blended with a weak nearby line at the blue edge, but 
we are not sure if this can account for the divergent abun- 
dance because the A604.60nm line seems to give a consis- 
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Table 1. Stellar basic parameters, kinematics, equivalent widths of the Si lines and abundances 
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tent abundance with that derived from the other S I lines 
for most of our sample stars. The loggf value presented in 
Kurucz (2002) (—0.78) will give an even higher abundance 
than that presented in Table ||. Therefore, we derived a 
differential loggf value for this line based on the spec- 
trum of HD693, which has a high signal-to- noise above 
300. We get loge Q (S) = 7.20 (excluding the A604.60nm 
line), which is the same as the meteorite value (Grevesse & 
Sauval 1998). We emphasize that the inclusion of this line 
or not in the abundance analysis for stars will not change 
the result of this work because it gives a consistent abun- 
dance with that derived from other lines. Experi menta l 
logg/ values for five Si I lines are taken from Garz ( 1973|) , 
and the differential loggf val ues fo r the remaining Si I lines 
are taken from Chen et al. (200C). The enhancement fac- 



tor, E*y, is assumed to be 2.5 for Fen and Si and 1.3 for 
Si I lines. There is no significant difference in the derived 
sulphur and iron abundances if is changed from 2.5 to 
1.5, because the Si and Fen lines are weak. Sil lines are 
somewhat stronger, but the maximum effect on [Si/Fe] for 
the most metal-rich star HD 1461 is only 0.05 dex when 
E 1 =2.h is adopted instead of 1.3. 

In sum, we obtained mean solar values: [Fe/H] = 0.02, 
[S/Fe] = -0.02 and [Si/Fe] = -0.01 relative to the me- 



teorite abundances given by Grevesse & Sauval ( 1998 ). 
Stellar abundances, presented in Table [l], are relative to 
these values, so that the solar ratios [X/Fe] in the following 
figures are exactly zero. 

3.3. Abundance uncertainties 

In order to investigate which lines, Fe I or Fe II, are the best 
indicator for the iron abundance in deriving the [S/Fe] 
ratio, we show in Table || the sensitivity of the derived 
abundances to the stellar parameters for a representative 
star HD 28620. This star has T cS = 6114 K, log g = 4.11, 
[Fe/H] = —0.38, which is typical for the range of temper- 
ature, gravity, and metallicity of our sample of stars. In 
these calculations, the logg f valu es for Fel lines are the 
same as those of Chen et al. ( 2000 ), and the deviations are 
obtained by a change of 100 K in effective temperature, 
0.1 dex in gravity, 0.1 dex in metallicity, and 0.3 kms -1 
in microturbulence. 

Clearly, sulphur abundances with respect to iron abun- 
dances derived from Fe n lines are more immune to vari- 
ations of atmospheric parameters as compared to sulphur 
relative to iron from Fel lines. This is primarily due to 
the high ionization potential of sulphur (10.36 eV), which 
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means that sulphur mainly occurs as S I atoms at temper- 
atures of —6000 K. Thus, we adopt the iron abundance 
derived from Fe II lines (hereafter the iron abundance) for 
all stars. Since the total error of [Si/Fe] ratio is only a bit 
larger when Si is measured relative to Fe derived from Fe II 
lines instead of Fel lines (see Table ||), both silicon and 
sulphur abundances are relative to iron abundances from 
Fen lines. A further advantage of using Fen lines is that 
non-LTE effects on [Fe/H] are bound to be small because 
nearly all iron occur as Fe n ions — a conclusion which is 
s uppo rted by detailed calculations by Thevenin & Idiart 
(|1999D . 

The line-to-line scatter of the deduced abundances 
gives an estimate of the uncertainty arising from the mea- 
sured equivalent widths. For silicon and iron, more than 
10 lines are available for most stars and the statistical er- 
ror of the mean abundance is less than 0.02 dex. With 2-5 
sulphur lines used in this work, the corresponding error of 
the mean abundance is of the order of 0.03 dex. Taking 
into account the errors from the uncertainties of the stel- 
lar parameters in Table ||, we estimate the total errors to 
be -0.06 dex for [Fe/H] and [Si/Fc], and -0.07 dex for 
[S/Fe]; the corresponding error bars are shown in Fig. 

The above error estimates do not, however, account for 
possible systematic errors introduced for example by the 
assumption of LTE or the use of ID hydrostatic model at- 
mospheres. The non-LTE effect on the S I A869.39nm and 
A869.46nm lines has been investigated by Takada-Hidai 
et al. (2002), who showed that this effect is negligible. For 
the other three S I lines, non-LTE effects are also expected 
to be insignificant because they have the same excitation 
potential as the A869.46nm line and they are weaker. As 
mentioned above, iron abundances derived from Fen lines 
are generally immune to non-LTE effects. Thus, non-LTE 
effects are expected to be negligible in the trend of [S /Fe] 
versus [Fe/H]. 

Following the same procedure as in Asplund et al. 
( |1999D and Asplund & Garcia Perez ( |200lD , we have esti- 
mated the effects of stellar granulation on the derived S 
and Si abundances using the new generation of ab-initio 
3D hydrodynamical model atmospheres (e.g. Asplund et 
al. 2000). A strictly differential comparison of the pre- 
dicted line strengths in the 3D models and ID MARCS 
model atmospheres with identical parameters has been 
performed for S I and Si I lines; the results for Fe II lines 
have been culled from Asplund & Garcia Perez (2001). 
As seen in Table [l| the ID results are relatively little af- 
fected due to the high excitation potential of the S I and 
Si I lines which forces them to be formed at similar atmo- 
spheric depths as the Fen lines. For the relevant metallic- 
ity range, the S/Fe and Si/Fe ratios differ by < 0.04 dex 
when using 3D model atmospheres compared to the ID 
case. To avoid interpolation within the still limited sample 
of granulation corrections and given their minor impact on 
the final results, we have made no attempt to apply such 
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Fig. 3. The [S/Fe] and [Si/Fe] ratio as a function of metal- 
licity. Filled circles denote the thin disk stars while open 
circles represent stars with kinematics typical for the thick 
disk. 

corrections to the ID abundances of the individual stars 
in the following discussion. 

4. Results and discussions 

4--1. Sulphur abundances vs. metallicity 

The [S/Fe] ratio as a function of metallicity is shown in 
Fig H It seems that the [S/Fe] ratio is nearly solar for 
[Fe/H] > —0.3, below which it increases with decreasing 
metallicity and an overabundance of —0.25 dex is indi- 
cated at [Fe/H] = —1.0. There may be a tendency for de- 
creasing [S /Fe] for super metal-rich stars but this is based 
on one star with [Fe/H] — +0.5 on ly. These re sults are 
in agreement with those of Francois ( 1987 , 1988 ), but the 



scatter in the [S/Fe] ratio at a given metallicity is signifi- 
cantly smaller in the present work. 



4-2. Reanalysis of data from Francois ( 198V , 198$ ) 



Israelian & Rebolo ( |2001 ) made a reanalysis of the data 
presented in Francois (1987, 1988) by using a new determi- 
nation of stellar parameters and suggest a linear increase 
of the [S/Fe] ratio with decreasing metallicity. One point in 
their paper is, however, strange: stars with [Fe/H] > —0.4 
have [S/Fe] around —0.3, which is significantly lower than 
the solar ratio. At the same time, they derived a solar 
sulphur abundance of 7.20, which is the same as the me- 
teoritic value. Furtherm ore, th e scatter in [S/Fe] obtained 
by Israelian & Rebolo ( 2001 ) is significantly larger than 
in the present work. 
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ments (open triangles) and our reanalysis of the Francois 
(p87L data (filled circles). 
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equivalent widths of Fe I lines published by Francois ( |1987| 



1988) were used to derive iron abundances. Adopted loggf 



values are -1.68 for Fei 867.48 nm and -0.38 for Fei 



869.95 nm taken from Nave et al. (1994); the very strong 
Fei 868.86 nm line was discarded because of its sensitivity 
to the uncertain damping constant. The solar equivalent 
widths for the Fe 1 867.48 nm and Fe 1 869.95 nm lines were 
measured from the solar flux spectrum by Kurucz et al. 
(1984). Again, stellar abundances (see Table ||) are relative 
to the Sun so that uncertainties from atomic line data can- 
cel out. The LTE iron abundances derived from Fei lines 
were corrected for non -LTE effects using the calculations 
of Thcvcnin & Idiart fll999| ). These non-LTE iron abun- 
dances are more consistent with those derived from Fe 11 
lines for our new sample of stars. As shown in Fig. |J, the 
scatter in the [S/Fe] vs. [Fe/H] trend at a given metallicity 
is quite small and the results from the reanalysis of the 
Francois data and the new measurements in the present 
work are consistent. 

Fig. H shows the deviations of [S/H] and [Fe/H] be- 
tween the present work and Israelian & Rebolo ( |2001 ) for 
14 stars in common. One of the 15 stars in our reanaly- 



sis sample was not included in Israelian & Rebolo (2001). 
It is clear that there is no significant deviation for the 
sulphur abundances between the two works, but the iron 
abundances show large differences . The iron abundances 
adopt e d by Israelian & Rebolo ( 2001 ) for 10 Frangois 
( 1987 , 1988 ) stars are significantly higher than the iron 
abundances in the present work, which explains the low 
[S/Fe] ratio (—0.3) they find at solar metallicity. We note 
that they adopt [Fc/H] values from different sources in the 
literature. Therefore, the iron abundances may be incon- 
sistent with T c ff and log g, and consequently also with the 
derived [S/H]. This explains both the large scatter and 
the offset with respect to the Sun in Fig. 3 of Israelian 
& Rebolo ( |200l| ), emphasizing the need for a proper and 
consistent treatment of the different species used to de- 
rive abundance ratios such as [S/Fe]. The adopted stellar 
parameters are, otherwise, quite similar for the two works. 



Fig. 5. Abundanc e com parison between this work and Is- 
raelian & Rebolo fl200lD for [S/H ] and [Fe/H] based on the 



reanalysis of the Frangois ( 1987 



stars. 



In order to investigate possible systematic differences 



relative to the work by Israelian & Rebolo (2001), we have 
made a reanalysis of sulphur and iron abundances for 15 
Frangois ( |1987| , |1988 ) stars with similar stellar parameters 
as our sample stars, i.e. 5500 < T e g < 6600 K, log g > 3.8 
and [Fe/H] > —1.1. The reanalysis was performed in the 
same way as that for our sample of stars, with parame- 
ters derived from Stromgren photometry and Hipparcos 
parallaxes. The microturbulence was esti mated from the 
relation presented in Edvardsson et al. ( |l993|) and the 



4-3. Correlation between sulphur and silicon abundances 

As described in the Introduction, one goal of this work is 
the comparison of sulphur with other a-capture elements 
in particular the closest element, silicon, being a good rep- 
resentative. The second reason for the choice of silicon is 
that we have many good silicon lines in the spectra and 
the oscillator strengths are quite reliable; abundances de- 
rived from different lines are consistent for both the Sun 
and stars. Thus we can expect that [Si/Fe] vs. [Fe/H] will 
be well defined. 

As shown in Fig. [| the trend of [Si/Fe] vs. [Fe/H] is 
similar to that of sulphur and the scatter at a given metal- 
licity is also the same. The one-to-one relation of [S/H] 
vs. [Si/H] is shown in Fig. 0; there is a very good cor- 
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Fig. 6. The correlation between [S/H] and [Si/H], and the 
[S/Si] ratio as a function of metallicity. 



relation between sulphur and silicon indicating common 
nucleosynthesis sites for these two elements. 

In order to investigate the dependence of the cor- 
relation of sulphur and silicon on metallicity, we plot 
the [S/Si] ratio as a function of [Fe/H] in the bottom of 
Fig. [| It seems that the sulphur abundance relative to 
silicon is constant for the whole [Fe/H] range investigated 
in this work. 



Goswami & Prantzos (2000) investigated the chem- 
ical evolution of many elements based on the yields of 



Woosley & Weaver (1995) for two cases: Case A consid- 
ered the constant (solar) metallicity yields while Case B 
corresponds to metallicity dependent yields. They found 
exactly the same trend for the sulphur and silicon evolu- 
tion in Case A and Case B in the whole metallicity range of 
-4.0 < [Fe/H] < +0.5. They suggested that the elements 
O, Si, S and Ca, behave as true primaries, whithout any 
metallicity dependence of their yields. The predicted ra- 
tios of S and Si in Goswami & Prantzos (2000 ) and other 
chemical evol ution models, e.g. Timmes et al. ( 1995 ), Chi- 
appini et al. ( 1997 ) and Samland et al. ( 1998 ), are consis- 
tent with the results of the present work for disk stars. 



4-4- Kinematics 

Stellar kinematics are useful information to interpret the 
stellar abundances and to probe the origin of the stars. 
We calculated the kinematics for our program stars based 
on proper motions and parallaxes in the Hipparcos survey 
and the radial velocity derived from the Doppler shift of 
the spectral lines in this work. The calculation of galactic 
space velocity (U, V, W) follows that presented in Johnson 
& Soderblom ( |1987 ). The correction of space velocity to 
the Local Standard of Rest is based on a solar motion, 
(-10.0, +5. 2,-1-7.2) kins" 1 in (U,V,W)as derived from 



Hipparcos data by Dehnen & Binney ( 1998 ). 

As shown in Table |l|, most of the sample stars are thin 
disk stars (Vlsr > -40kms _1 ). Exceptions are HD 10453 
and HD 22879, which have thick disk kinematics. The 
[S/Fe] ratio is about 0.2 dex for both stars, which is higher 
than most other stars. These data seem to be in line 
with the result by Edvardsson et al. (1993), who found 
an abundance gradient of a elements at the metallicity of 
[Fc/H] ~ —0.7. But the metallicities of the two thick disk 
stars are also the lowest, and thin disk kinematics does 
not imply a low [S/Fe] ratio. For example, HD 77967 and 
HD 159307 have thin disk kinematics, but both [S/Fe] and 
[Si/Fe] are close to the enhanced value of ~0.2 dex. There- 
fore, it seems that these abundances ratios are primarily 
a function of metallicity and the kinematics plays a minor 
role. 



5. Conclusions 

We have carried out new measurements of sulphur and sil- 
icon abundances for 26 main sequence stars in the metal- 
licity range of —1.0 < [Fe/H] < +0.5 and have established 
the trends of [Si/Fe] and [S/Fe] with a small scatter at a 
given metallicity. The result indicates a very good correla- 
tion between sulphur and silicon and the correlation does 
not depend on the metallicity. The results provide support 
for the chemical evolution models of Goswami & Prantzos 
( |2000D , Timmes et al . Q1995D , Chiappini et al. ( |1997| ) and 
Samland et al. ( 1998 ) for these two elements in the case of 
disk stars. It is concluded that the nucleosynthesis of sul- 
phur is a proxy for silicon, which is the nearest a-capture 
element of sulphur in the periodic table of elements. 
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Table 2. The atomic line data for Si, Sil and Fell lines, 
and solar flux equivalent widths and abundances. Note 
that the mean solar sulphur abundance is derived exclud- 
ing the A604.60nm line. 



Ele 


A 


Xi 


log gf 


ref. 




EW 


Abu. 




nm 


eV 








mA 




Si 


604.603 


7.87 


-0.51 


1 


2.5 


14.5 


7.32: 


Si 


605.267 


7.87 


-0.63 


2 


2.5 


9.8 


7.25 


Si 


675.717 


7.87 


-0.31 


2 


2.5 


16.1 


7.21 


Si 


869.396 


7.87 


-0.52 


2 


2.5 


11.0 


7.19 


Si 


869.464 


7.87 


0.05 


2 


2.5 


28.9 


7.17 




Mean loge Q (S) = 


7.20 ± 0.03 






Si I 


566.556 


4.92 


— 2.04 


3 


1.3 


40.9 


7.57 


Sil 


E C! f\ C\ AO 

569.043 


4.93 


—1.87 


3 


1.3 


53.8 


7.62 


Si i 


570.111 


4.93 


—2.05 


3 


1.3 


37.1 


7.52 


Si i 


570.840 


4.95 


—1.40 


1 


1.3 


76.0 


7.49 


Sil 


577.215 


5.08 


— 1.66 


1 


1.3 


54.8 


7.55 


Sil 


579.308 


4.93 


—1.95 


1 


1.3 


47.2 


7.59 


Sil 


579.786 


4.95 


—2.05 


3 


1.3 


43.1 


7.64 


Sil 


594.855 


5.08 


— 1.19 


1 


1.3 


87.1 


7.57 


Sil 


612.503 


5.61 


— 1.54 


1 


1.3 


33.3 


7.54 


Sil 


n-\ A O/lfl 

014.^49 


5.62 


—1.48 


1 


1.3 


38.7 


7.59 


Sii 


ol4.5Uz 


5.62 


— 1.43 


1 


1.3 


39.6 


7.55 


Sii 


703.491 


5.87 


—0.81 


1 


1.3 


69.6 


7.57 


Sii 


TOO <?01 


5.61 


—1.30 


1 


1.3 


47.6 


7.53 


Sii 


740.579 


5.61 


—0.68 


1 


1.3 


93.6 


7.54 


Sii 


741.596 


5.61 


—0.71 


3 


1.3 


91.5 


7.54 


Sii 


791.838 


5.95 


-0.54 


1 


1.3 


95.3 


7.65 


Sii 


793.235 


5.96 


-0.35 


1 


1.3 


113.6 


7.64 


Sii 


872.802 


6.18 


-0.36 


1 


1.3 


91.2 


7.53 






Mean lo; 


je (Si)= 


=7.57+0.03 






Fen 


599.138 


3.15 


-3.61 


4 


2.5 


31.9 


7.52 


Fe ii 


608.410 


3.20 


-3.86 


4 


2.5 


20.7 


7.52 


Fe ii 


611.332 


3.22 


-4.16 


1 


2.5 


11.3 


7.50 


Fen 


614.925 


3.89 


-2.77 


4 


2.5 


37.5 


7.50 


TT 

pen 


617.940 


5.57 


-2.65 


4 


2.5 


o. z 




Fen 


623.995 


3.89 


-3.49 


4 


2.5 


13.3 


7.53 


Fen 


624.756 


3.89 


-2.38 


4 


2.5 


56.4 


7.52 


Fen 


638.371 


5.55 


-2.32 


4 


2.5 


7.2 


7.54 


Fen 


641.693 


3.89 


-2.79 


4 


2.5 


40.7 


7.59 


Fen 


643.268 


2.89 


-3.63 


5 


2.5 


41.4 


7.49 


Fen 


645.639 


3.90 


-2.13 


4 


2.5 


65.2 


7.46 


Fen 


651.608 


2.89 


-3.38 


1 


2.5 


54.0 


7.50 


Fen 


722.240 


3.89 


-3.35 


4 


2.5 


18.9 


7.57 


Fen 


722.446 


3.89 


-3.29 


4 


2.5 


19.7 


7.54 


Fen 


744.933 


3.89 


-3.22 


4 


2.5 


21.1 


7.51 


Fen 


747.970 


3.89 


-3.64 


4 


2.5 


9.7 


7.49 


Fen 


751.583 


3.90 


-3.48 


4 


2.5 


13.1 


7.50 


Fen 


771.173 


3.90 


-2.59 


4 


2.5 


49.3 


7.54 



Mean loge Q (Fe)=7.52±0.04 



Table 3. Dependences of relative abundances on the stel- 
lar parameters. 







HD28620 






T eff =6114 K, log g= 


=4.11, [Fe/H] = 


=-0.38, 6 


=1.6 kms" 1 




T c ff 


log g 


[Fc/H] 






+100 K 


+0.1 dex 


+0.1 dex 


+0.3 km s" 1 


[Fei/H] 


0.07 


0.00 


0.00 


-0.05 


[Feii/H] 


-0.01 


0.04 


0.01 


-0.04 


[S i/Fei] 


-0.12 


0.03 


0.00 


0.05 


[S i/Feii] 


-0.04 


-0.01 


-0.01 


0.04 


[Sii/Fei] 


-0.04 


0.00 


0.00 


-0.03 


[Sil/Fen] 


0.04 


-0.04 


-0.01 


-0.02 



Table 4. Comparison of the Si, S and Fe LTE abun- 
dances derived with ID hydrostatic and 3D hydrodynam- 
ical model atmospheres. The 3D abundances are those 
which reproduce the equivalent widths computed using 
a ID model atmosphere with identical parameters, a mi- 
croturbulence of fturb = 1.0 km s^ 1 and the abundances 
given in the third and fourth columns 



T e a log g 


[Fc/H] 


a/Fef 


A (S/Fe) 


b A (Si/Fe) c 


[K] [cgs] 










5767 4.44 


+0.0 


+0.00 


-0.02 


-0.01 


5822 4.44 


-1.0 


+0.40 


-0.03 


-0.04 


6191 4.04 


+0.0 


+0.00 


+0.02 


-0.01 


6180 4.04 


-1.0 


+0.40 


+0.01 


-0.03 


a The a-enhancement 


for the 


S and 


Si abundances. 


For these 


calculations 


, loge(S) 


=7.33, lo 


ge(Si)=7.55 and 


log e(Fe) = 


7.50 have been assumed for the Sun following 



Grevesse & Sauval (1998). 
b Based on the five Si lines listed in Table 2 and the six 

Fell lines used by Asplund & Garcia Perez (2001). 
c Based on the four Sil 570.1, 594.8, 722.6 and 872.8 nm 

lines and the six Fe II lines used by Asplund & Garcia Perez 

(2001). 



References for log g f v alues: (1) differe ntial values in this 
work ; (2) Kurucz ^002| ; (3 ) Gar z |l973| ; (4) Biemont et al. 
199l] ; (5) Hannaford et al. |l992| 
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Table 5. Stellar parameters and relative abundances for 
the rcanalysis of data from Frangois (1987, 1988). 



HD 



K 



log 



6 

kms 



[Fe/H] 
LTE 



[Fe/H] 
NLTE 



[S/H] 
LTE 



09984 


0896 


o no 
3.93 


t.() 


— 0.88 


— 0. (4 


—0.0/ 


63077 


5825 


4.15 


1.2 


—0.89 


—0.75 


—0.57 


69897 


6227 


4.20 


1.6 


-0.50 


-0.41 


-0.41 


76932 


5870 


4.08 


1.3 


-0.98 


-0.83 


-0.63 


88218 


5661 


3.94 


1.4 


-0.53 


-0.44 


-0.26 


91324 


6123 


3.95 


1.8 


-0.60 


-0.50 


-0.49 


102365 


5562 


4.39 


1.0 


-0.39 


-0.32 


-0.35 


106516 


6135 


4.30 


1.5 


-0.79 


-0.66 


-0.55 


136352 


5584 


4.27 


1.0 


-0.58 


-0.48 


-0.44 


139211 


6231 


4.12 


1.7 


-0.26 


-0.21 


-0.16 


148816 


5784 


4.07 


1.4 


-1.03 


-0.87 


-0.54 


157089 


5712 


4.00 


1.4 


-0.79 


-0.66 


-0.44 


193901 


5662 


4.48 


1.0 


-1.13 


-0.96 


-0.67 


201891 


5825 


4.21 


1.3 


-1.21 


-1.03 


-0.86 


203608 


6094 


4.29 


1.4 


-0.91 


-0.77 


-0.69 



